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^\ Abstract 
O ■ 

' In an external constant magnetic field, so strong that the electron Larmour length is much shorter 

C/3 ! 

than its Compton length, we consider the modification of the Coulomb potential of a point charge 

^r~i . owing to the vacuum polarization. We establish a short-range component of the static interaction 
> ■ 

' in the Larmour scale, expressed as a Yukawa- like law, and reveal the corresponding "photon mass" 

' parameter. The electrostatic force regains its long-range character in the Compton scale: the tail 

rn ; 

. of the potential follows an anisotropic Coulomb law, decreasing away from the charge slower along 

o : 

the magnetic field and faster across. In the infinite-magnetic-field limit the potential is confined 

o ■ 

' to an infinitely thin string passing though the charge parallel to the external field. This is the 

^ ■ first evidence for dimensional reduction in the photon sector of quantum electrodynamics. The 

b : 

5t , one-dimensional form of the potential on the string is derived that includes a (5-function centered 

in the charge. The nonrelativistic ground-state energy of a hydrogenlike atom is found with its use 
and shown not to be infinite in the infinite-field limit, contrary to what was commonly accepted 
before, when the vacuum polarization had been ignored. These results may be useful for studying 
properties of matter at the surface of extremely magnetized neutron stars. 
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I. INTRODUCTION 



The fact that the vacuum, in which an external magnetic field is present, is an opti- 
cally anisotropic medium, has been known, perhaps, since the time when the nonlinearity of 
quantum electrodynamics was first recognized: in a nonlinear theory electromagnetic fields 
do interact with one another, provided that the strength of at least one of them is of the 
order of or larger than the characteristic value Bq = m^/e ~ 4.4 x lO^^G, where m and 
e are electron mass and charge, respectively. [Henceforth, we set /?. = c = 1 and refer to 
the Heaviside-Lorentz system of units.] If the external field is strong, other fields interact 
with it, the result of the interaction depending upon the direction specified by the external 
field, hence the anisotropy. Depending on the wave amplitude, the electromagnetic waves 
propagation in this medium may be considered as a nonlinear process [l], including the 
transformation of one photon into two [2j or more photons, or taken in the linear approxi- 
mation with respect to the amplitude. In the latter case, the second-rank polarization tensor 
is responsible for the properties of the medium. In the kinematic domain where the photon 
absorption processes like electron-positron pair creation are not allowed, the polarization 
tensor is symmetric and real, and the medium is transparent and birefringent [sj. In the 
absorption domain the medium is dichroic {4 1 . The limit of low frequency and momentum 
belongs to the transparency domain and corresponds to a constant anisotropic dielectric 
permeability of the medium. In this limit the polarization tensor may be obtained by dif- 
ferentiations with respect to the fields of an effective Lagrangian, calculated on the class 

n 

of constant external electric and magnetic fields. For small values of these fields [5[ and 
for extremely large {4] fields the polarization operator was in this way considered using the 



effective Lagrangian of Heisenberg-Euler calculated 



a] within the one-loop approximation. 



(The two-loop calculations are also available The knowledge of this limit is useful for 
studying the dielectric screening of the fields that are (almost) static and (almost) constant 
in space. For more general purposes, however, this limit is not sufficient, and one should cal- 
culate the polarization tensor directly, using the Feynman diagram technique of the Furry 
picture in the external magnetic field. On the photon mass shell, i.e. when the photon 
energy, ko, and 3-momentum, k, are related by the free vacuum dispersion law = k^, 
such calculations were done by Adler and Constantinescu i9i|. The results obtained are 
appropriate for handling the photon propagation in weakly dispersive medium, when the dis- 
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persion law does not essentially deviate from its vacuum shape. The polarization operator 
for the case of general relation between the photon mass and momentum was calculated by 
Batalin and Shabad [l^, Tsai 11], Baier et al. [3], and Melrose and Stoneham 13]. This 
gave the possibility of studying the photon propagation [l4] under the conditions where the 
deviation from the vacuum dispersion law may be very strong either due to the phenomenon 
of the cyclotron resonance in the vacuum polarizatio n |l5(] - this phenomenon is responsible 



for the effect of photon capture by a magnetic field [l6 



much larger than Bq 



19 



21 



or due to magnetic fields, 



|22| , or due to the both circumstances (see Ref . [23| where 
the photon capture effect was extended to low frequencies for extra-large fields). 

Although much work has been devoted to study of electromagnetic wave propagation in 



the magnetized vacuum, problems of e^ 
attract sufficient attention, save Refs. 



ectro- and magneto-statics in this medium did not 



24 



251 ]. where corrections to the Coulomb law were 



found when these are small: for B/Bq -C 1 in 2J], or at large distances from the source 
for 1 ^ B / Bq ^ Sna^^ in 25|, where a = e'^/An = 1/137 in the fine-structure constant. 
Here we proceed with an investigation of some electrostatics in the presence of a strong 
external magnetic field in the vacuum to find that for sufficiently large b = B/Bq ^ I the 
electric field produced by a pointlike charge at rest may be significantly modified by the 
vacuum polarization, the modification being determined by the characteristic factor ab. We 
note first of all that expressions for the dielectric permeability of the magnetized vacuum 
obtained from the Heisenberg-Euler Lagrangian are applicable only as far as the fields slowly 
varying in the space are concerned. Otherwise, the spatial dispersion becomes important. 
For this reason, when considering the electric field produced by a pointlike electric charge 
in the present paper, we address again to the polarization tensor calculated off-shell in 
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12 



13 



141 ]. taking it in the static limit /cq = 0, but keeping the dependence on 



k. The corresponding spatial dispersion effects will be essential for getting some important 
features of the modified Coulomb potential. 

Using the tensor decomposition of the polarization operator and the photon Green func- 



tion, established in 



10 



14| in an approximation-independent way, we find that photons of 



only one polarization mode (mode-2 in the nomenclature of these references, see below) may 
be carriers of electrostatic force. This is in agreement with the fact that the electromagnetic 
field of these photons is, in the static limit, purely electric and longitudinal. The photons of 
the other two modes mediate in this limit the magneto-static field of constant currents. 
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In magnetic fields B ^ Bq, which we are deahng with when describing the static field, 
produced in the magnetized vacuum by a point electric charge, the electron Larmour length 
Lb = (c-B)^^/^ = Ac&^^^^ is much less than the electron Compton length Ac = 
Therefore, two different scales occur in the problem: the Larmour scale and the Compton 
scale. 

A simplifying expression for the mode-2 eigenvalue of the polarization operator is used, 
valid for such fields. It was first obtained by Loskutov and Skobele y | 25 | within a special 
two-dimensional technique intended for large fields, and by Shabad 19], and Melrose and 
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I2I, Q, Q in 



Stoneham 13] as the asymptote of the mode-2 eigenvalue calculated [10 
the one-loop approximation (see 23(] for the detailed derivation of the large-field asymptotic 
behavior.) The most important, now widely accepted, fact about this asymptotic behavior 



(see, e.g., the monographs 



20 



26 



27( 1) is that the mode-2 eigenvalue contains a term linearly 
growing with the magnetic field, seen already jj] if one deals with nondispersive small mo- 
mentum approximation, inferrable from the Heisenberg-Euler Lagrangian. It is sometimes 
expected that this term - it appears in the denominator of the photon propagator and hence 
of the expression for the potential - should lead to suppression of the interaction mediated 
by mode-2 photons. In a different problem we already had an opportunity to show that this 
is not exactly the case 28| . In 23| the impact of the linear term on the refractive index was 
considered. In the present paper we demonstrate that this term is also crucial for the most 
important features of the potential of a point charge. 

Correspondingly to the two scales inherent in the problem, the potential is divided into 
two additive parts, out of which the first one, called short-range, decreases exponentially at 
distances r of a few Larmour lengths from the source, but retains the usual q/Airr singularity 
near the origin r = 0, where the charge q is placed. (The anisotropy shows itself no sooner 
than in the third term of the Laurent expansion around the singular point r = - see 
Appendix I.) The second part, called long-range, slowly decreases away from the charge 
following an anisotropic Coulomb law, but remains finite close to the charge. The linear 
term mentioned above, is responsible for the fact that a scaling regime of the short-range 
part occurs, characterized by a comparatively simple universal function, independent of 
the magnetic field. This function gives the potential of a point charge in the energy units 
of Lg^ = (eBy^'^ as a dimensionless function of the space coordinates in the units of the 
electron Larmour length Lb. Excluding the closest vicinity of the charge, its form coincides 
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with the Yukawa law [see Eq. (125|) below] characterized by the dimensionless mass parameter 
2a I n (which is the topological mass of the two-dimensional Schwinger electrodynamics 29|, 
measured in inverse Larmour lengths). Thus, this mass governs the exponential (isotropic) 
decrease of the potential away from its source at the distances, large in the Larmour scale. 
In the formal limit of infinite magnetic field the short-range part becomes the 5-function 
with its center in the charge. As one moves farther from the source, the Yukawa decrease 
ceases, and the potential coincides with its long-range part. It approaches, for distances 
from the charge, large in the Compton scale, the anisotropic Coulomb shape of the form 
of Eq. ( l38l) (that might have been also derived disregarding the spatial dispersion). The 
law of decreasing along the magnetic field is unaffected by the latter, the decrease is the 
fastest in the direction orthogonal to the magnetic field. Thus, the linearly growing term 
in the mode-2 eigenvalue of the polarization operator leads, first, to the faster decrease of 
the potential in the direction across the magnetic field for large distances, and, second, to 
its steeper shape for small distances. This may be recognized as suppression, indeed. On 
the other hand, the long-distance behavior along the magnetic field, as well as the standard 
Coulomb singularity near the source 30|| do not sense the magnetic field at all, no matter 
how strong it is. 

Perhaps, the most interesting feature of the potential produced by a point charge is that, 
as the external magnetic field tends to infinity, the whole potential becomes concentrated 
inside an infinitely thin string that includes the charge and is directed parallel to the magnetic 
field. The electric lines of force produced by the charge are gathered inside the string. 
The string is the b = oo limit of an ellipsoidal equipotential surface. The potential along 
the string as a function of the longitudinal distance from the charge is just the infinite- 
magnetic-field limit of the long-range part of the potential (see Fig. [6] in Section IV) plus 
the (5-function contribution from the short-range part. The string potential first grows 
with distance logarithmically and linearly (starting with negative values) in-between the 
Larmour and Compton distances and hence provides "confinement" in this scale. For the 
string formation, again, the above-discussed term, linearly growing with the magnetic field, 
is responsible. To conclude about its presence, a consideration of the Heisenberg-Euler 
Lagrangian might have been sufficient. However, for calculating the string potential, the 
effect of spatial dispersion was important. In contrast to the quark-antiquark string in 
the lattice QCD, the string potential of the present paper stops growing after reaching the 



6 



Compton distances from the charge and approaches zero following the Coulomb law l/Anr 
in accord with the fact that the infrared singularity in QED is milder than in QCD and 
insufficient to provide the infra-red custody. 

The appearance of the string notifies the reduction of QED to two dimensions (one time- 
one space) in the photon sector in the infinite-magnetic- field limit, which implies a new result 



3l| . because what was known before was the reduction to two dimensions in the electron 



sector. The latter circumstance is a common knowledge and is well understood referring 
to the fact that electrons are confined to the lowest Landau level, so only one degree of 
freedom - that along the magnetic field - survives to remain dynamical. For instance. 



it was demonstrated in 



281], that the Bethe-Salpeter equation describing the interaction 



between electrons and positrons acquires in this limit fully Lorentz-covariant form in the 
two-dimensional space with the metrics (1,-1). 

Analogously, the nonrelativistic Schrodinger equation for an electron in the field of the 
nucleus of a hydrogenlike atom is known to become a differential equation with respect 
to the longitudinal distance between the two particles. Unless the vacuum polarization is 
taken into account, the standard result 32] is that due to the singularity of the Coulomb 
potential in the origin, the ground-state energy in this problem tends to negative infinity as 
the magnetic field unlimitedly grows. The conclusion of the present paper is that if the string 
potential obtained is used in the Schrodinger equation the ground state energy remains finite 
just because the singularity of the string potential in the origin has the 5-function character. 

Another conclusion concerns the critical nucleus charge Zcr making the threshold of its 
instability manifested in spontaneous free positron production. The known fact here 33] 
is that becomes reasonably smaller than = 137 for large magnetic fields. This 
result depends on the same unboundedness from below of the energy spectrum of the Dirac 
operator caused by the same Coulomb singularity of the static potential. Therefore, it should 
be revised if the vacuum polarization of the Coulomn potential is taken into account. 

The paper is organized as follows. In Section II we present an approximation-independent 
form of the potential of a point charge in terms of the mode-2 component of the photon 
propagator and define the division of the potential corresponding to the one-loop polarization 
operator in the asymptotical region of large magnetic fields into the short- and long-range 
parts. In Section III we consider the short-range part as determined by an expansion near the 
universal function corresponding to the scaling regime. It is obtained as a one-fold integral. 
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We also establish the (5-function character of the short-range part in the limit of infinite 
magnetic field. In Section IV the anisotropic Coulomb law is obtained for large distances for 
the long-range part of the potential by applying mathematical means, differnt for different 
remote regions in the space. The limiting, b = oo, form of the long-range part is studied for 
short and long distances on the string. In Section V we estimate the limiting value of the 
ground state energy of a hydrogenlike atom by considering the Schrodinger equation with the 
vacuum-polarization-modified potential and using the shallow-well-potential method. Also 
a perturbation correction to the ground state valid for the fields in the range 1 <^h <^2Tx/a 
is found. A role the radiative modification of Coulomb potential may play when the Dirac 
equation is used is discussed. In Section VI we briefly discuss possible applications of our 
results to physics of strongly magnetized neutron stars. In Appendix I serving Section II 
we derive the asymptotic expansion of the potential near the point r = and study the 
coefficients in this expansion as functions of the magnetic field. Also an analog of Uehling- 
Silber correction to the Coulomb potential valid in the interval 1 ^ 6 ^ 27r/Q; is derived. 
In Appendix II we deal with a simplified potential that models the short-range Yukawa-like 
potential and also has the 5-function as a limiting form. In this case explicit solution of 



3^. The 



the Schrodinger equation can be obtained with the use of the method of Refs. 
finiteness of the ground energy is demonstrated. 

Some results of the present paper were shortly reported in our previous paper [35|. See 
also the preliminary publication 36|], more detailed in certain points. 



II. GENERAL REPRESENTATION FOR THE STATIC POTENTIAL OF A 
POINTLIKE CHARGE: ONE-LOOP APPROXIMATION IN THE LARGE MAG- 
NETIC FIELD DOMAIN 

Electromagnetic 4- vector potential produced by the 4-current ju{y) is 

A,{x) = J D^,{x - y)ny)d''y, /i, z/ = 0, 1, 2, 3. (1) 

Here x and y are 4-coordinates, and D^y{x — y) is the photon Green function in a magnetic 
field in the coordinate representation. The metric in the Minkowski space is defined so that 
diag g^i, = (1,-1,— 1,-1). Eq. ([1]) defines the 4- vector potential with the arbitrariness of a 
free-field solution, including the gauge arbitrariness. If the photon Green function is chosen 
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as causal, only the gauge arbitrariness remains. 

The current of a pointlike static charge q, placed in the point y = is 

fiy)=qS.oS'iy). (2) 

Hence 

/oo 
Df,o{xo -?/o,x)dyo 
-OO 

roc poo 

= q Df,o{xo + yo,:)i.)dyo = q D^o{yo,^)dyo- (3) 



This 4-vector potential is also static. 

If there is no magnetic field, and the photon propagator is free and taken in the Feynman 
gauge (with the pole handled in a causal way) 

D,Ax -y) = D%{x -y)^ ■,^J{J_yY 
only the zeroth component of the 4-vector potential is present: 

This is the Coulomb potential in the Heaviside-Lorentz system of units used throughout. 

Let there be an external magnetic field B directed along axis 3 in the Lorentz frame 
where the charge q is at rest in the origin x = , and no external electric field exists in this 
frame. Call this frame special. Define the Fourier transform as 

Df,t,{x) = j exp{ikx)D^^{k) d'^k, /x, z/ = 0, 1, 2, 3. (6) 

Then ([3]) becomes 

^/.(x) = J exp[i{koyo - kx)]Df,oik) d'^kdyo = J D^o(0, k) exp(-ikx) d^/c. (7) 



The four 4-eigenvectors a = 1,2, 3, 4, of the polarization tensor U^^, [lO|, [l^, ll5| take 



in the special frame (and arbitrary normalization) the form - the components are counted 
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downwards as z/ = 0, 1, 2, 3 



b(2) 



/ \ f ko\ 

+ kl 
k2 k2 

\oJ \k,J 

( \ 

k2 
-ki 

\ y 

Among them there is only one, whose zeroth component survives the substitution k^ = 0. 
It is bt/ . This imphes that out of the four ingredients of the general decomposition of the 
photon propagator 



b(4) 



fkA 




\ko/ 

ki 
k2 



(8) 



D^,{k) = Y,'^aik) 



a=l 



(bW)2 ' 



Va{k) 



~[k^ + K,{k)]-\ a={l,2,3} 
arbitrary, a=4 , 



(9) 



where Ka{k) are scalar eigenvalues of the polarization tensor: 

u;{k)\>l^^ = ^,{k)\>^^\ K,{k) = o., 



(10) 



only the term with a = 2, P2(/s)b|?'*b[/^'*/(^*"^'')^5 participates in (JTj), i.e. only mode-2 (virtual) 
photons may be a carrier of electro-static interaction, and not photons of modes 1 and 2, 
nor the purely gauge mode 4. Bearing in mind that (b*^^^)^ = ~ ^O' have 



An fx) 



-ikx 



d^k 



27r)3 y k2 - K2i0,klk 



2 7.2 \' 



Al,2,3(x) = . 



(11) 



Here k\ = k\ + k^. Thus, the static charge gives rise to electric field only, as it might be 
expected. The gauge arbitrariness in the choice of the photon propagator Di{k)\) 



(4), (4) 



T)/^{}i)k^ky indicated in ([9]) has no effect in (ITT]) . Certainly, the potential (ITT]) is defined up 
to gauge transformations. 

The result that only mode-2 photons mediate electrostatic interaction may be understood, 
if we inspect electric and magnetic components of the fields of the eigen-modes obtained from 
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their 4- vector potentials ([8]) in the standard way: e*^") = koh^"''^ — kbg'^'' , h^"^ = k x b^"^ . These 
are Q 

e« = -^fco, hW = X k3), (12) 

k± k± 



ej^ = kxfca, eg 



f = |(*i-*2), hf^> = -*o (kx X I) , (13) 



where the cross stands for the vector product, and the boldfaced letters with subscripts 3 and 
-L denote vectors along the directions, parallel and perpendicular to the external magnetic 
field, respectively. 

The photon energy and momenta here are not, generally, related by any dispersion law. 
Therefore, we may discuss polarizations of virtual, off-shell photons - carriers of the inter- 
action - basing on Eqs. f|T2|) - f|T4|) . The electric field e in mode 1 is parallel to k^, in mode 
2 it lies in the plane containing the vectors k, B, in mode 3 it is orthogonal to this plane, 
i.e. mode 3 is always transversely polarized, e*^^^k = 0. For the special case of the virtual 
photon propagation transverse to the external magnetic field, k^ = 0, (this reduces to the 
general case of propagation under any angle ^ 7^ by a Lorentz boost along the external 
magnetic field), mode 2 is transversely polarized , e^^^k = 0, as is always the case for mode 3. 
Mode 1 for transverse propagation, k^ = 0, is longitudinally polarized, e*-^-* x k = 0, and its 
magnetic field is zero. The lowest-lying cyclotron resonance of the vacuum polarization jisl , 
the one that corresponds to the threshold /cq — fc| = 4m^ of creation of the pair of electron 
and positron in the lowest Landau state each, belon gs t o mode 2. It gives rise to the photon 



capture effect with the photon turning into a free [16|] or bound 17|, [iSj electron-positron 
pair. Another consequence of the cyclotron resonance is that a real photon of mode 2 un- 
dergoes the strongest refraction in the large magnetic field limit even if its frequency is 
far beyond the pair production threshold. 

In the static limit ko = the magnetic field in mode 2 disappears, h^^^ = 0, while its 
electric field is collinear with k, e^^^^ = k. It becomes a purely longitudinal virtual photon. 
Unlike mode 2, in modes 1 and 3 in the static limit ko = only the magnetic fields survive: 
e(i,2) ^ = X B, h^^^ = -k^ks, hf^ = k\, h(i'3)k = 0. (Here normalizations 
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are arbitrary and kept fixed only within the same mode). A virtual mode-1 photon carries 
magneto-static interaction. It is responsible for magnetic field produced by a current fiowing 
through a straight-linear conductor oriented along the external magnetic field. In accordance 
with the above formula for h^^^ its magnetic field is orthogonal both to B and to the radial 
direction in the transverse plane k^, along which the magnetic field of the current decreases. 
The mode-3 photon contributes as an interaction carrier in the problem of a magneto-static 
field produced by a solenoid with its axis along axis 3. In the present paper, however, we 
do not consider magneto-static problems. 

In the asymptotic limit of high magnetic field eB ^ A;|, B ^ rr? je = Bn the eigenvalue 



/t2(0,k), as calculated within the one-loop approximation |10l . I 111 . 1 121 . Il3l . Il4| . with the 



Q, Q, Q,J3, 



accuracy of terms that grow with B only as its logarithm and slower is 



19| 



.,o,.i.i)^-^e.p(-^)r(A), 

where b = B/Bq and 

Note that < T{y) < 1 for ?/ e (0, oo). It will be demonstrated in the subsequent sections 
that the asymptote T{y — > 0) ~ 2?//3 in ( |T6l) is responsible for the large-distance Coulomb- 
like behavior of the potential in the direction orthogonal to the field, while the asymptotic 
value T(oo) = 1 introduces a sort of photon mass and governs the short-range Yukawa-like 
part of the potential [see Eqs. ([24]), ([25]) and ([38]) below]. 

Other eigenvalues, ki^s, do not contain the coefficient b that provides the linear increase 
of K2 ([T5]) with the magnetic field. Therefore, in the polarization tensor, whose covariant 
decomposition is 

n,.(A:) = 5^/..(A:)-^, (17) 

a=l ^ ' 

the components fi, u = (0, 3) dominate in the high magnetic field limit in accord with the 
idea about the two-dimensioning of the photon sector. 

Expression for k,2 ([151) is to be used in ([TT]) or, equivalently, in the expression 

q f°° ( /"°° e~''^^^^dA;3 \ 



12 



N -2 



<< -3 




FIG. 1: Electron energy — e^o(2;3,0) in the modified Coulomb potential (llSp of the charge q = Ze 
plotted along the axis X3 passing through the charge q parallel to the magnetic field, = 0. Thin 
solid lines correspond to four values of the magnetic field (from left to right): b = 10^, b = 10^, 
b = 10^ and b = W^. Bold solid line is the Coulomb law j4q (X3, 0) = q/{4Trx3). Thin lines approach 
asymptotically the bold line at the both edges of the figure. Thick dashed broken line corresponds 
to 6 = 00. The abscissa represents the distance in the units (Im)"^. The ordinate represents the 
potential in the units 2Zam = Z x 7.46 keV. 



obtained from (fTTj) by integration over the angle between the 2- vectors and x^, which are 
projections of k and x onto the plane transverse to the magnetic field. In f|T8l) x± = \/ x\ 
and Jo is the Bessel function of order zero. We explained in 35|] why the fcs-integration here 
may be extended up to the value j/csl = cx) in spite of the limitation on the validity of (ITSll 
indicated above. 

The results of computer calculation following Eq. (fT8|) of the shapes of the potential for 
the two cases, x_l = and X3 = 0, are given in Figs. [T]and[2l respectively js^. The curves in 
Figs. [Hand [2] manifest the standard Coulomb singularity l/|x|, when |x| 0. Then, they fall 
up rather sharply, following a Yukawa-like law within the range of several Larmour lengths 
Lb to reach the asymptotic long-range regime that is the Coulomb law y4Q(x3,0) = g/47ra;3 
for = and |x3| ^ mT^ and Ao(0,x^) = g/47rx'j_, x'^ = x±^yl + ab/37r, for X3 = and 
Xj_ ^ [see Eq. f l38|) below]. In what follows we shall be commenting on the features of 
the computed curves referring to analytical considerations. 

The nontrivial - other than the leading asymptote ~ kl near the point = k± = - de- 



pendence of the polarization operator eigenvalue (1151) on the photon momentum components 
ks, k± is the spatial dispersion. We shall see in Section 3 that it is important in the vicinity 
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0.2 0.4 0.6 0.8 1 1.2 

x,[ (2m)-i] 



FIG. 2: Electron energy — eAo(0,xj_) in the modified Coulomb potential (jlSp plotted along the 
axis x± passing through the charge q transversely to the magnetic field, X3 = 0. Conventions are 
the same as in Fig. [TJ Thin lines approach asymptotically the solid line at the left lower edge of 
the figure and the short dotted lines (f38]) Ao(0,x±) = q/{4:-Kx'j_) at the upper right edge. 

of the charge, where the field has a large gradient. As for the anisotropic behavior far from 
the charge, to be studied in Section 4, only the above asymptote is essential, inferable also 
from the Heisenberg-Euler Lagrangian. 

In Appendix 1 we consider the singular asymptotic behavior of the potential in the 
vicinity of the point charge and present its expansion near |x| = 0. Now we shall consider 
separately two additive parts of the potential that decrease by different ways with increase 
of the distance x. 

It is useful to subdivide identically the potential f|T8l) as: 



Ao(x) = As.r.(x) + Ai.r.(x) 



(19) 



with 



and 



^.r.(x) 



Ar.(x) 



2(27r) 



Jo(A;_lx_l) 



A;2 + H + 



exp 



X 



/c^ ~l~ /C3 ~l~ 



exp 



7^ ( ^3 ) 



p + p ^ 2a6rni 



fc2 



dfc 



exp 



2m^b 



Eq. (!20l) is the potential (fTT!) or (fTSj) taken with the substitution T 
inside K2 (fT5l) . We shall see in what follows that r.(x) is a Yukawa-like potential, singular in 



(20) 



dkl. (21) 



T(oo) 
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the origin, that exponentially decreases at distances of about ^Ji: j^a ^ 15Lb, while y4i.i..(x) 
is a finite function that slowly decreases at large distances (greater than the Compton length 
m~^) following what will be called anisotropic Coulomb law. This is the reason why we shall 
call Eq. fl20p the short-range part and Eq. fl21l) the long-range part of the potential. 

Consider first the short-range part (120|) and the limiting form it takes in the infinite-field 
limit. 



III. SHORT-RANGE PART 



A. The scaling regime 



It is remarkable to note that the short-range part of the potential fl20l) . as measured 
in the inverse Larmour length Lg^ = V eB units is a universal, magnetic-field-independent 
function of coordinates measured in Larmour units Lb- To establish this scaling regime let 
us make the change of variables in the integral fl20l) = k^L^, k± = k±L-B and define the 
new dimensionless coordinates X3 = x^L-q, x± = x^L-q. Then Eq. ( 120|) becomes 

^-.(x) = ttttA^ / Uk^x^) / ^ r.x - (22) 



2(27r)2LB 7o i-00 kl + kl + ^ exp (- 



This is an even function of X3. The fca-integration here can be performed by calculating the 
residues in the poles on the imaginary axis in the points 



±i 



A^i + -exp(-4l (23) 



TT \ 2 

with the upper sign taken for X3 > and the lower one for X3 < 0. Finally one gets 



As.r(x) = ^ = ^ / Jo(fcA)^^= dfcj. (24) 

Here the universal function A(x) depends on the magnetic field through its spatial arguments 
x only. Eq. fl20|) (or fl22l) and (|2^ ) is illustrated in Fig. [3] drawn for a;^ = by a computer. 
The simple representation (|2^ can be further simplified if X3 or x± are large in the Larmour 
scale: l^sl ^ 1, or \x_[] ^1. In this case the integration in (1241) is restricted to the domain 
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0.2 0.4 0.6 0.8 1 1.2 
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FIG. 3: Electron energy —eAs.rX^s^^) i^i the short-range part Eq. ([20l) of the potential plotted 
against the (absolute value of) longitudinal distance for x_l = (thin dashed lines from left to 
right correspond to b = 10®, 10^, 10^, 10^). The thick dashed broken line corresponds to 6 = oo. 
All the rest is the same as in Figs [1] and [2j 



k'j_ <^ 1 where the exponential exp(— ^5^/2) should be taken as unity. Then flMl) is reduced 
to the Yukawa law for the short-range part of the potential 



A ~ /lYr^^ _ 9 exp[-(2a/7r)VV^+x^] _ q exp [-(2a&/7r)^/^m|x| ] 

47rLB ^yxl_ + xj 47r |x| 



It reflects the Debye screening of the charge by the polarized vacuum. Eq. (1251) can be 
established if we return to the previous representation (122|) . which can then be traced back 
to (flHl) with 



«2(0,cx),0) = ^T(oo) = m2— (26) 



substituted for — ^2(0, A;^) in the denominator. Here M is the "effective photon mass" 



of Ref. 



38|. Write it as 



M = . (27) 

71V2 Lb 

The deviation of (12^ from (l25l) when x± and X3 are both small is not very important against 
the background of the singularity of the short-range part of the potential y4s,r.(x) near the 
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charge provided by the divergency of (120!) in the origin x=0. (See Appendix 1 for the 
asymptotic expansion of the potential near the origin). Therefore, the Yukawa law fl25l) is 
approximately fulfilled also in the close vicinity of the charge. Correspondingly, the curves 
in Fig. [3] could not be distinguished from the Yukawa law fl25l) in the scale of the figure. 

Eqs. fl24p and (123]) establish the short-range character of the static electromagnetic forces 
in the Larmour scale. The corresponding effective mass (!26|) coincides with the topological 
photon mass esch/ in the 2-dimensional Schwinger electrodynamics 29| provided that the 



dimensional fermion charge esch of that theory is identified as esch = \/7r/2e/LB. Stress that 
the zero photon mass understood as its rest energy is also present as a consequence of the 
gauge invariance reflected in the approximation- independent relation ^^(O, 0, 0) = 0. Corre- 
spondingly, the potential, produced by a static charge, should be long-range for sufficiently 
large distances. This is the case, indeed. The long-range character of the electromagnetic 
interaction is restored at the distances of the Compton scale, as we shall see in the next 
Section. The carrier of the long-ranged interaction will be fl2Tl) . The Debye screening 

obtained here in the vacuum completely depends on the fact that the function ( fT6l) tends to 
unity for large longitudinal momentum fcs oo, i.e. on spatial dispersion. In this point the 
situation is different from the case of a medium, where the Debye screening can be achieved 



391] in expressions, obtained from the thermodynamical potential, which is the analog of 



the effective Lagrangian for that case. The difference with the medium is also in that the 
long-range part of the static potential is absent in that case in spite of the gauge invariance, 
because it implies that the appropriate polarization tensor components should disappear in 
the long-wave limit ko = 0, k=0 only if k is set equal to zero first [4^, thus providing the 
zero value to the photon magnetic mass. 



B. The limiting b = 00 form 

The short-range part fl20l) and fl22l) tends to zero, when b ^ 00 for any nonzero distance, 
X 7^ 0, from the charge. This follows from the fact that A{x), defined in (12^ . tends to 
zero with the exponential speed fl25l) when X3 = x^/L-q or x±_ = x±/Lb tends to infinity. 
As the magnetic field b grows more and more, the curves representing the potential (l24ll 
for the special case of = in Fig. [3] stick closer and closer to the vertical axis, the 
spacing between the curves and this axis becoming infinitely thin in the limit b = 00. The 
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area {q/2Ti)S of the region restricted by any curve (12^ and the xs-axis in the domain 
— oo < X3 < —Lb, Lb < 0:3 < 00 



•2 -L2eiexp 



^5 = 2/ As.r.(x3,0)dX3 = 2 / A(J3,0)dJ3 = f" / "3 r^T^ (28) 



2 

is a finite number, with 5* = 2.180, independent of the magnetic field. If the Yukawa 
approximation (125|1 is used in (128!) in place of (12^ . approximately the same numerical value 
for S is achieved: 5 = —'E^.{ — ^J2a/'K) ~ 2.176, where 

/oo 
exp(-y)y-My (29) 

is the exponential integral. In the limit 6 = oo. Lb = the width of the strip |x3| < Lb 
excluded from the integration in Eq. (l28l) is zero, and the latter becomes the whole area 
above the limiting potential. Thus, in the infinite-magnetic-field limit the short-range part 
f l20|) of the potential taken on the axis drawn through the point charge q along the magnetic 
field direction becomes the 5-function: 

^.r.(x3,0)|,=^ = 2.180:f5(x3). (30) 

The limiting 5-function here is understood in the following sense. Given a test-function 
t(x3), one has 

-Lb 



lim 

»oo 



As.r.(a;3,0)t(x3)dx3 + / As.r.(a;3,0)t(x3)dx3 

Jlb 



2.180^t(0). (31) 
27r 



[This equation directly follows from the scaling law, the first equality in (Elj)]. In this sense 
it will be used in Section IV and Appendix 2, where we shall see that the 5-singularity of 
the potential (130|) leads only to a finite contribution to the atomic ground-state energy in 
an infinite magnetic field in contrast to the contribution of the primary Coulomb potential. 

Analogously, we may write a 5-function for the limiting form of the short-range part of 
the potential along any direction x, |x| = in the plane orthogonal to the magnetic field 
containing the point charge q. In place of ( |28l) one has 



—S^ = 2 [ A, r (0, x)dx = 2 [ A{0, x)dx 
2vr Jl^ Ji 



'Lb 

q_ r 1 - fc±Lo(fc±) + 4^[Lo(fc±)Hi(fc^) - Ji(fc^)Ho(fc^)] ^ q_^ ^^^^ 
'-J^ /^i + ^exp(-|) 
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Here Jo,i and Ho,i are, resp., the Bessel and Struve functions of orders zero ^and one. We used 
the integral 6.512.8 and the representation 8.551.1 in the reference book [41,] for calculating 
. Note that the Struve functions at large argument decrease and oscillate asymptotically 
like the Neuman functions; besides Hi includes a constant asymptotic term 2/7r. The integral 
f l32p converges: the divergence caused by the unity term in the nominator is cancelled by 
the two products of two oscillating asymptotic terms in the square brackets. 

From fl32l) and the fact that for 7^ the short-range potential flMl) disappears in the 
b ^ 00 limit we have 

A.r.{0,x)\,^^ = 2.l78^S{x). (33) 



The different coefficients in f l5U]) and ([3SD manifest the anisotropy. Note that the Cou. 
singularity of the (short-range part of) the potential in the origin g/47r|x| is isotropic 



omb 



30|. 



IV. LONG-RANGE PART 



We have finished the consideration of the short-range part and will proceed with consid- 
ering the long-range part Air (x) fl2Tl) . Simplifying expressions will be obtained for large- 
distance behavior in Subsection A and for the long-range part taken on the axis x_l = in 
Subsection B. We shall see in Subsection B that in the limit b = 00 the long-range part, as 
well as the whole potential, is concentrated on this axis, making an infinitely thin tube or 
string. We shall study the potential along the string in more detail in Subsection B. 



A. Long-distance behavior of the long-range part ^i r.(x) 

Once we have seen in the previous Section that the short-range part yls.r.(x) is as a matter 
of fact concentrated within the region of a few Lb, for larger distances, |x| > m~^, the whole 
potential flTSl) and its long-range part are the same. For this reason in this Subsection 
we shall deal directly with (ITSIl . 



1. Large x± in Larmour scale 



For large transverse distances the term linearly growing with the magnetic field (fT5!) leads 
to suppression of the static potential in the transverse direction. 
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To be more precise, consider the region 

Lb /„.s 

Once the Bessel function Jq in (ITSl) oscillates and decreases for large values of its argument 
k±Xj_, the main contribution into the integral over k\ in f|T8|) comes from the integration 
variable domain <^ 2m?h, and the dependence upon k^ in ^2 niay thus be disregarded. 
Then the ^^."iiitegration in ( ITSl) can be explicitly performed to give (we use Eq.6. 532.4 of 



the reference book 
^o(a^3,a: 



m 

2q 



^0 kl - ^2(0, A;i,0)^ cos(A;3X3)d/c3, (35) 



- (27r)2 

where /Cq is the McDonald function of order zero, and 

K,(0,40) = -H^m^T(A). (36) 

As the McDonald function /Co decreases exponentially when its argument increases, only 
small values of the square root contribute into integral (135|1 . and the /cs-integration domain 
in it is restricted to the interval /c| <^ 4m^, wherein 

T ( ~ (37) 

V 4m^ / 6m? 



Then the potential form ( l35ll becomes (we use Eq. 6.671.14 of the reference book 4l| ) 



^0(2:3, x±) ~ [ /Co ( x±k3\ 1 + ^1 COs(fc3X3)dfc3 



(27r)2 7o - ■ Svr 

= X_L I 1 + — , > X_L. (38) 



Eq. ( l38il is an anisotropic Coulomb law, according to which the attraction force de- 
creases with distance from the source along the transverse direction faster than along 
the magnetic field (remind that h = (B/Bq) ^ 1), but remains long-range. The 
equipotential surface is an ellipsoid stretched along the magnetic field. The electric field 
of the charge E = —VAo{xs,Xj_), as written in Cartesian components, is the vector 
(g/27r)(xi + /32a;i)-3/2 {/S^xi (3'^X2 X3), where /? = (1 + ab/Srry/^ It is not directed 
towards the charge, but makes an angle (p with the radius- vector r, cos(f) = {x^ + l3'^x\){x'^ + 
/3^a;5^)^^/^(x| + x^)~^/^. If x_l 7^ 0, in the limit of infinite magnetic field, P 00, the electric 
field of the point charge is directed normally to the axis 0:3, since the ratio {E^/E^) 

20 



(although Es and E± are both equal to zero in this limit outside the string). But if x± = 0, 
the electric field is directed along the external magnetic field. It looks like the electric field 
compresses the string. This regime corresponds to the dielectric permeability of the vacuum 
independent of the frequency, with its dependence on k (spatial dispersion) being reduced 
solely to that upon the angle in the space (c/, |23l|). 

The result (l38l) is in agreement with the curves in Fig. [2] in the large x± domain. 



2. Large 

It remains to consider the remote coordinate region of large 0:3, complementary to ([3i 
To this end we apply the residue method to the inner integral over in (fTSll . Using 
the integral representation f|T6|) the function k,2 f|T5|) may be, for a fixed positive value of 
kj_, analytically continued from the real values of the variable into the whole complex 
plane of it, cut along two fragments of the imaginary axis. In the lower half-plane the cut 
runs from Im k^ = —2m down to Im k^ = —00, while in the upper half-plane it extends 
within the limits 2m < Im A;3 < 00. Other singularities of the A;3-integrand in (fT8l) are 
poles yielded by zeros of the denominator, i.e. solutions of the equation (associated with 
the photon dispersion equation) 

kl + kj- K2{0,klkl) = 0. (39) 



As k± varies within the limits (0, 00) two roots of this equation fc^ = ±iK(k±) move a. 



ong 



the imaginary axis from the point K(0) = to the points kf = ±iK(oo) = ±i2m jl4l.ll5l.l23|. 
There is yet another branch of the solution to equation fl39p . corresponding to the photon 
absorption via the 7 — e"'"e~-decay, but the corresponding poles lie in the nonphysical sheet 
of the described complex plane, behind the cuts, and will not be of importance for the 
consideration below. 

Let us consider positive values of X3. Negative values can be handled in an analogous way. 
Turning the positive part of integration path < A;3 < 00 clockwise to the lower half-plane 
by the angle tt/2, and the negative part —00 < k^ < counterclockwise by the same angle, 
and referring to the fact that the exponential exp(— 1/^3X3) decreases, for X3 > 0, in the lower 
half-plane of k^ as {k^l ^ 00 so that the integrals over the remote arcs may be omitted, we 
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get a representation for the inner integral in (fTS!) 



'-00 kl + kl-^,{0,klkl) 

/"OO 

. / e-^'"'^'''A{\k3f,kl)d\ks\-i2TTexp[-K{kl)x3]Res{kl), (40) 

J2m 



where Res^kjJ designates the residue of the expression D2{0,—\ks\'^,k'j_) = 
[kl-lksl"^ - K2{0,-\kl\,kl))~^ in the pole k^ = -iK{kl), while A{\ks\^kl) = 
1^2(0, — l^sp + iO, fc^) — 1^2(0, — l^sP — iO, k'^) is the cut discontinuity. It was explained above 
that < KlkjJ < 2m everywhere but in the limit k± — » 00, where K = 2m. Consequently 
the residue term in fHOl) dominates over the cut-discontinuity term everywhere in the 
/c^-integration domain in the outer integral in (fTSll . except for the region near k± = 00. 
In this limit, however, k,2 disappears due to the exponential in ( fTSi) . together with the cut 
discontinuity, since the latter is only due to the branching points in the function (fTGll . 
Therefore, keeping the residue term in fBUl) as the leading one, we neglect the contribution 
that decreases with large longitudinal distance at least as fast as exp(— 2m|x3|). In this 
way we come to the asymptotic representation of the potential (fT8|) in the region of large 
longitudinal distances jxsl (2m) ~^ (negative values of 0:3 at this step are also included - 
to handle them one should rotate the fragments of the integration path in the directions 
opposite to the above) 

, ^, g r '^o(^^^^) i-K{ki)\x3\] dkj 

""'^""'-SnJ, Kikl)[l + Hi-K^ikl),kl)] ' ^^'^ 

where 

../,9 ,9x 2abm^ ( k^\ \ d ^ ( k\\ 

ml *i) = — -p (-^) d^i^ (i;^ j ■ (42) 

iere K'^i}?^ is the solution of equation fl5^ in the negative region of the variable k\ - see 



231] for its form. K[oo) = 2m, K{0) = 0. T is given by (fT6l) . 

Due to the exponential factor in the integrand of (14T!) . for large X3 the main contribution 
comes from the integration region of k± that provides minimum to the function K{k±). The 
minimum value of K{k±) is zero. It is achieved in the point A;^ = - a manifestation of the 
fact that the photon mass defined as its rest energy is strictly equal to zero owing to the gauge 
invariance: Ka(^o = A;3 = k = 0) = 0. In view of fl36p and (13 7p . near the point k± = the 



dispersion equation ( l39l) has the solution K{k±) = k±/^yi + ab/37T. Simultaneously, near 
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the minimum point 1 + H{0, 0) = l + ab/Sn. With these substitutions and the use of 6.611.1 



of 



4l[ |. Eq. pri) becomes again the anisotropic Coulomb law fl38l) (q/An) /[{x'j_Y + a;?!^/^ 



We have, therefore, established its validity everywhere in the region remote from the center, 
irrespectively of the direction. 

In agreement with this result the curves in Fig. [T] for Aq{x3,0) approach the Coulomb 
law q/4:n\x3\ as |x3| grows. The difference between the potential Ao{x3,0) and its large-xa 




FIG. 4: Electron energy —eAQ{x3, x±) in the modified Coulomb potential (j4ip with q = Ze 
plotted against the transverse coordinate x± with the longitudinal coordinate fixed at the large 
value X3 = 10(2m)~^. Thin solid line corresponds to the magnetic field value B = W'^Bq. Bold 
solid line is the standard Coulomb law ([5]) — e^Q (xs, x±) = —2aZm[{2mx ±)'^ + 100]~^/^. The thin 
line is indistinguishable from the anisotropic Coulomb law (j38p in the scale of the drawing. The 
coordinate axes are the same as in Fig. [2] 



asymptote q/4n\x3\ decreases in Fig. [T] at least as fast as exp(— 2m|x3|) (see 36| for the 
derivation of the latter statement). 

Eq. fHTj) was used for computer calculation with the large value 0:3 = lOm^^. It has 
led to the curve shown in Fig. |H In the region (l34l) it agrees with the result (1381) . valid in 
that region [Lb = 0.02(2m)^^ for b = 10^]. In practice fj4T]) and fl38|) are the same. A small 
deviation of the potential curve Ao{10/2m, x±) from (1381) may be seen in Fig. 5 of Ref. j36| . 
drawn in a more detailed scale for small x±. 
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B. The long-range part on the axis x_l = and its limiting form for b = oo 

Curves drawn for Aij-X^s, 0) by a computer following Eq. (12T!) are presented in Fig. O 
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FIG. 5: Electron energy —eAiT- (xs,0) in the long-range part ()2ip of the potential for b = oo 
(dashed thick line) and b = 10®, 10^, 3 x 10'^, 10^ (dashed lines from bottom to top) 



Here we study the form the long-range part (!2T]) of the potential takes in the limit b = oo. 

First consider the case x± 7^ 0, X3 7^ 0. As we saw in Subsection B of Section 3 the 
short-range part of the potential tends in this case to zero as 6 ^ 00. Therefore, the limits 
of the whole potential and of its long-range part are the same. For this reason to achieve the 
claimed goal it is sufficient to consider the limit of (|T8|l . After the change of the integration 
variable kj_ = k^m\/h Eq. (ITSl) becomes 



An fx) 



2(27r) 



Jo{k±m\/hx±) 



e-ifc3X3(ifc3dfc2 



I + ^ + — exp 



(43) 



When b ^ {Ti/2a) one can disregard the ratio in the denominator. 

For any finite x±_ the argument of the Bessel function in f H3l) is large, therefore we may 
use the same procedure as the one that led us from f|T8l) to fl35|) and fl38l) . Then we obtain 

2q I . fob 



Ao{X'i,XA_)\ 



6— >oo, x^^O 



i2ny 



/Co x±k 




cos{ksXs)dk3 



qm 



47r ^x\{ah/?>Tx)+xl A^JaTx/?>xi_ 



(44) 
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This means that outside the xs-axis the potential (fT8|) turns to zero as the ratio Lb/xj,. 
Since its short-range part fl20|) or fl25l) decreases with b exponentially, the result (jH]) holds 
for the long-range part (121 p as well. 

The situation is different on the axis x± = 0. By making in Eq. fl^Tl) the same change 
of the variable k± as above and, again, neglecting k^/Am'^b in the denominators we come to 
the limiting (6 = oo) form of the long-range part of the potential x± = 0, independent of 
the magnetic field 



^ roc r<yo 



^0 



1 



^i + ^exp(-|)T(^) A;i + ^exp(-| 



dkl. (45) 



This is the analytic representation of the envelope curve in Fig. [H To understand this, 
note that the overall potential is the sum of the short- and long-range parts, according to 
( IT9|) . Therefore by combining the curves in Figs. [3] and [5] we come to the pattern presented 
in Fig. [6l which is the detailing of Fig. [H Each potential curve drawn for a certain value 
of the magnetic field approaches, as the distance from the charge along the xs-axis grows, 
the corresponding (dashed) curve transferred from Fig. O But even prior to this, the latter 
approaches the thick dashed curve, which is the common envelope of the curves in Fig. [5] 
and the whole potential curves in Figs. [1] and [61 

We continue by studying the long-range part of the potential along the string, Eq. fHSl) . 
To separate the part independent of the fine-structure constant a the internal integral here 
is integrated by parts to yield: 

d fci + ^exp(4i)r(y)~ 

In ■ 



dki 



dkl kl+^fexpl^^^) /l + l 

kl + ^exp(^)T{y) 
+ -r^ln ^)iy . (46) 

Then fj45l) becomes 



2io ^i + ^y fci + ^exp^z^ 



2a ( 

q p cos{k3X3)dk3dkl 1^ K ^ ~ V~ 



T (ii) - 1 



2(2-PA Jo f,^~^y \ A;i + ^exp(4i 
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(47) 
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FIG. 6: Four electron energy — eAo(x3, 0) curves in the modified potential (jlSp for b = 10^, 10^, 3 x 
10^, and 10^ (thin solid lines from left to right) approaching their corresponding long-range parts 
— e^i.r.(2;3, 0), Eq. (pT|) . shown in Fig. [5] (four dashed lines from bottom to top). Thick dashed 
broken line corresponds to the string potential Aq(x3,0) I^^q^. Its vertical fragment symbolizes the 
5- function ([3T]) . 



where the first term (it is worth recalling here that within the integration limits T{y) is a 
positive function, lesser than unity) 

1 



Ar.(a;3,0)| 



6=00,0=0 



cos{ksXs) In 



is independent of a, whereas the second term behaves as 



T 



47r3 



-a In a) 



cos(A;3X3) 



T 



4^2 



^3 
4^2 



dfe. 



(48) 



(49) 



when a tends to zero, i.e. is nonanalytic in a = 0. The reason for the nonanalyticity 
and for the nonvanishing of Ai r. (2^3, 0)1 6=00 0=0 that that a chain of diagrams has been 
as a matter of fact summed when solving the Dyson- Schwinger equation that led to the 
expression for the photon Green function with the one-loop polarization operator (fT^ 
substituted into it. In the result flTT]) thus obtained the two limits b = 00 and a = are not 
permutable. 

Eq. (HHI) may be referred to as a fitting approximation for the envelope (H5|l . simpler than 
( H7|) . It is presumably useful for making rough estimates with the accuracy to {—a/n) In a = 
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0.011. It might have been obtained if the exponential exp(— in 0451] had been merely 
replaced by unity. The integral fHSj) is converging at the both integration limits due to the 
asymptotic properties of the function T{y) indicated below its definition (fT6l) and represents 
a function that decreases at large longitudinal distances following the Coulomb law {q/Airxs). 

The limiting curve (H5ll . or P7|) . (148!) for the long-range part of the potential fl2T]) crosses 
the axis X3 = in the point Ar.(0, 0)|f,^^ = (1.4240 - 0.0088) (gm/27r) = 1.4152(gm/27r). 
Here two numerical contributions from the first and the second terms in (1471) are presented 
separately. It is intriguing how close the numerical coefficient in the expression for the 
intercept of the envelope and the = axis is to V2 = 1.4142. More precise value of 
would be achieved by the infinite-magnetic-field limit of the long-range part of the potential 
in the point where its charge is located, if the fine-structure constant were 1/121. Higher-loop 
calculations may improve this figure. 

Identifying the above-calculated a-dependent coefficient 1.4152 supposedly with a/2, an 
interesting observation would follow: if the charge q is taken equal to the electron charge 
e, {Z = 1), the string potential undergoes the increment between the point X3 where the 
charge is located and the infinitely remote point = 00, equal to 

AAu.{x3) ^ A.r.(0,0) 1,=^ - Ai.,.(oo,0) 1,^^ = , (50) 

nyZ Ac 

where Ac = Tn~^ is the Compton length. This "work function over a unit charge" differs 
from the photon mass (|27|) in that the Larmour dimensioning has been replaced by the 
Compton one. 

In the interval jxal <^ (2m) the envelope curve fj45l) looks roughly in Fig. [T]as a linearly 
growing potential, the same as this is believed to be the case for the confining potential along 
quark-antiquark string in QCD in the limit of zero lattice spacing. We may say that in QED 
the "confinement" occurs within distances smaller than the Compton length, whereas for 
larger distances - thanks to the fact that the infrared behavior in QED is weaker than in 
QCD - the growth of the potential ceases and it approaches the zero value along the Coulomb 
asymptote. As a matter of fact the growth of the potential is only nearly linear. 

To estabhsh its true character consider the difference Ai j-X^s, O)!^^^^ — A.r.(0, 0)1^,^^ and 
change to the new integration variable u = k^x^ in the integrals ( HHl) and ( 1471) . Then the 
argument of the function T{y) becomes {u"^ /Am'^xl) and should be considered as large, once 
Arn^xl <^ 1. According to Eq. (1T6|) for large y one has [T(?/) — 1] ~ (— ln2y/2?/). As long as 
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this tends to zero with y oo, we may substitute ln(l + x) ~ x, x <^ 1 for the logarithms 
in (HHj) and fH7|) . In this way we obtain 



where 

^-^n j7:2 

1 



ln(2m|x3|) - -ln2 + 7 - 1 



,(51) 



/(«) 



exp ( ] dfc^ 



1 + 



fc2 



.2 I 2q 



-exp 



/ 



137.036 



4.533, /(a)l^^o ~ -Ina, 



and 7 = 0.577 is the Euler constant. We have made use of the two standard integrals 



(COSM — 1) , TT 

^ dw = 



'0 ^ JO 

Finally, for small distances along the string, 2|x3|m <^ 1, the long-range potential has the 
form: 



2' 



^ - Inw dw = — (7 — 1). 



(53) 



A.r.(a;3,0)| 



6=00 



qm 
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[1.4152 + 0.495 X 2m|x3|(ln(2m|x3|) - 0.77)], 



(54) 



This should be additively combined with the 5-function (1301) . to which the short-range part 
^s.r.{^3,)^) is reduced in the limit h = 00, to form the string potential. It is this potential 
that is responsible for forming the spectrum of an atom in the strong magnetic field, to 
consideration of which we are proceeding. 



V. RADIATIVE SHIFT OF ELECTRON GROUND-STATE ENERGY IN A HY- 
DROGENLIKE ATOM IN A STRONG MAGNETIC FIELD 

In this section we shall study how the ground-state energy of a hydrogenlike atom at rest 
in a strong magnetic field is modified by the radiative corrections to the Coulomb potential 
considered above. 

The nonrelativistic electron in an atom, whose nucleus has the charge q = Ze,is described 



by the one- dimensional Schrodinger equation 32] 



1 d^^ixs) 
2m dx| 



- eAo{x3, x± = 0)^(X3) = E^{x3), {x^l > Lb = {eB)-'/\ (55) 



if the atom does not move transverse to the magnetic field - which is the case as long as we 
are interested in its ground state. The one- dimensional Schrodinger equation (j55|) is valid 
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in the region |x3| > Lb and is efficient provided that Lb ^ (Ib/Z, where ob = 1/ma is the 
Bohr radius. 

If the unmodified Coulomb potential (jl]) taken at = is used for Aq{xs, 0) in equation 
f lSS]) {q = Ze henceforward), the ground-state energy value 

Eo = -2ZVmln2^, (56) 



that follows 
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34| from equation fl55|) is unbounded from below, in other words tends 
to negative infinity as the magnetic field grows. The reason is that the one-dimensional 
Coulomb potential is too singular, the singularity being regularized by the Larmour radius. 
In equation (1551) this regularization acts as the cut-off of the definition region |x3| > L-q of 
equation (!55l) . The regularization is lifted by letting h tend to infinity. Lb 0, and hence 
Eq —>■ — oo. On the contrary, the radiation corrections studied here, yielded the conclusion 
that the Coulomb g/(47rx3) singularity of the one-dimensional potential in X3 = had 
been changed to the (5-function ( !30l) . This sort of singularity is not expected to cause 
an unboundedness of the energy spectrum. To confirm this, we solve in Appendix II the 
Schrodinger equation flS^ with a potential that models the short-range part As ^..(^s, 0) fl2U]) 
alone and also tends to 5-function as 6 ^ 00. The resulting ground-state energy approaches 
in this limit the finite, magnetic-field-independent value given by Eq. (11051) of Appendix II. 
The genuine level must be significantly lower due to impact of the long-range part of the 
potential Air,{x3,0) (j2Tl) shown in Fig. [5l 

A. Extremely large magnetic fields 

To estimate the ground-state energy E in the limiting case 5 — »• 00 we apply here 



the shallow- well approximation of Ref. 42], appropriate since the electron potential 



V = —eAo{x3,x± = 0) has a small depth {\V\ ^ (ma^) ^, where a is the range of the 
forces in the well). In this case, the value of E may be estimated as 

E^-2m(^j^ eAo(x3,0)cix3^ . (57) 

Here it is taken into account that the electron potential is symmetrical, Aq{—x^,Q) = 
Ao{x3,0). 
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To find first the contribution of tlie long-range part into (157|) rewrite (12T|) as (g = eZ, 

k±m\/b = k±, 2mk = ) 



2 roc poo 

eAij.{x3,0) = 2aZm— / / cos(2m/cx3) 



dA; 



kl + lk-^ + 'f exp ( ) T (A;2) kl + f F + ^ exp ( 



dfc 



(58) 



Tlie b = oo limit of this expression is Eq. or fl47l) . In the problem under consideration the 
potential falls following the Coulomb law, and hence, according to [42], the upper integration 
limit in (1571) should be replaced by the Bohr radius ob = {ma)~^. Then the contribution of 
the long-rang part fl58l) into the ground-state energy is determined by the integral 



,1 



oo POO f2k\ 1 
ev4i.r.(a;3,0)dx3 = aZ- I / sin — - 
Jo Jo \a J k 

1 



1 



/ci + |P + ^exp(-^ 



dkdkl. (59) 



Analogously, the contribution of the short-range part fl20l) into the ground state energy 
according to (!28l) is the magnetic-field-independent constant 



eAs,,,(a;3,0)da;3 = aZ 2.18. 



As a matter of fact, only the contribution of the first, a-independent term 



(60) 



eAr. (3^3,0) 



b=oo, Q=0 



dxi 



-azi /""sinf—VnT(A;2)^ = aZ 6.392 (61) 



'0 J ' ' k 

is important, whereas the second term in (H7|l only corrects the value 6.392 in the third 
decimal number (+6 ■ 10^'^). Combining fl6Tl) with fl28l) we get from fl57|) the finite value of 
the energy level of a hydrogenlike atom in the limit b = oo 



Ey, 



lim 



-2mZ'^a^ 73.6 



-Z^ X 4.0 keV. 



(62) 



This result reproduces with good accuracy the value obtained by us earlier [35'] with the use 
of a graphically fitted formula in place of Eq. ( 1471) . 

The Loudon-Elliott energy (!56|) would overrun the limiting energy (162|) already for the 
magnetic field as large as 6 = 6600, when yet the proton size R ~ lO^^^cm remains much 
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smaller than the Larmour length, R <^ Lb- The ground level reaches 92% of its limiting 
value for 6 = 5 x 10^. After the magnetic field reaches the value b = 1.5 x 10^, when R and 
Lb equalize, the Coulomb potential is cut off at the proton size, X3 = R. Setting L-q = R 
in (1561) we would get the minimum value for the Loudon-Elliott energy {Z = 1) to be —5.6 
keV, which is essentially lower than (1^2]) . 



B. Moderate magnetic fields 

For moderate magnetic fields lying in the range 1^6^ 10^ the additive radiative 
correction to the Coulomb law, as calculated in Subsection B of Appendix I keeping the first 
power of ab/2'K in the power series expansion of Ao(x3, 0), 

A^ofe.O) ^ »^ Texp i-"^) cos'* (63) 
Stt^ Jo V COS0 y 

may be considered as perturbation. Therefore, the radiative shift to the ground-state energy 
evel can in this case be calculated using the purely Coulomb (normalized) wave function 
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^>{x,) = exp (-P^) (64) 



as unperturbed. Here 60 is the "quantum defect" for the Coulomb problem 

^o = ilnf-|^l (65) 



Z yAa^Z"^^ 

Calculating the average of fl63l) multiplied by e with the wave function (j6lj) we find the 
perturbation to the Loudon-Elliott ground-state energy {q = Ze ) 

° ~ 27raB5o Jo m + ^ ~ 2n X cos0 + V«' 

One sees that for magnetic fields within the scope of applicability of the expansion in powers 
of ab, 1^6^ lO'^, where (!63l) is valid, the quantum defect 60 ^ a, or (l/Za) — 8.454 ^ 
In^b/Z"^), provided that Z < 11. Then, for the light hydrogenlike atoms Eq. (|66|) can be 
further simplified to 

E-Eo = In ^ = Z% 0.18 (^ In A + 1) eV. (67) 

" Svr Aa^Z^ ^8.454 Z^ ' ^ ^ 
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C. Relativistic corrections to electron motion 



The value (162|) makes about 1% and more of the electron rest mass, hence the question 
about relativistic corrections may arise. 

When the Dirac equation with stationary Coulomb potential is considered in infinitely 
growing magnetic field, the effect of unlimited lowering of the energy level down to — cx3 
is enhanced as compared to the Schrodinger equation due to the known fact [5] that the 
potential is squared after the Dirac equation is reduced to one-component second-order 
differential equation. Therefore, we should face an one-dimensional second-order equation 
with the stronger singularity (aZ/xa)^, apart from the singularity aZ/x^ already present in 
f lSS]) . For this reason one may expect that the ground-state energy would tend to negative 
infinity faster than the logarithm squared in fl56|) . Anyway, according to the (numerical 



part of the) analysis in Ref. [33[, it rather sharply approaches the border of the lower 
continuum E = —m, where the instability with respect to free positron production opens, 
analogous to what happens (without any magnetic field) for nuclei with the supercritical 
charge Z > 137 (to be more precise, Z > 170, once the finite size of the nucleus is taken into 
account) 43j. Whereas for infinite magnetic field the unlimited sinking of the level occurs 
already for infinitesimal Coulomb attraction Za — 0, for large, but finite magnetic field the 
Coulomb-induced quadratic singularity is cut off at the Larmour length, hence the ground 
level reaches the lower continuum at finite Za. The dependence of the corresponding critical 
value of Z on the magnetic field was found long ago by Oraevskii, Rez and Semikoz [33|, who 
claimed, for instance, that already for the values of magnetic fields b = 10^ to 10^ (that may 
exist near neutron stars according to the estimates available at present time, see below), the 
critical value of the nuclear charge lies within the reasonable range in the Periodic Table 
Z = 55 to Z = 90. This result is to be reconsidered now that we have established the 
important alteration in the singular behavior of the modified Coulomb potential proved to 
be crucial for the Schrodinger equation with huge magnetic field. 

The next level of relativistic description of atomic (or positronium) spectrum based on 
the static potential would be that via the Bethe-Salpeter equation with the so-called equal- 
time Anzatz wherein the recoil of the point source of the electrostatic field (a nucleus or 
a positron) is taken into account, but the retardation effects in the relative motion of the 
electron and the nucleus (positron) are disregarded. The corresponding results established 
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m 



should be also subjected to revision. (This statement does not concern the 



conclusions about the effect of photon capture through positronium formation in the pulsar 
magnetospheres made in 

The matters stand differently when very deep relativistic effects are dealt with. The latter 
come into play for magnetic fields tens of orders of magnitude higher than those for which the 
asymptotic limit in the present context is saturated (i.e., than, say, b = 10^°). Retardation 
effects make the static potential an insufficient quantity to take on the responsibility for 
forming bound states, since the full electromagnetic interaction is mediated by all the three 
photon modes in (Q. Unlike fT5l), the polarization operator eigenvalues ni^s of two other 



modes do not include 



13 



19 



20 



25( 1 the fast-growing factor b, and hence the interaction 



singular on the light cone — = characteristic of the free photon propagator is not 
suppressed in these modes. Correspondingly, the infinite deepening of the energy level. 



considered in our papers [28| for positronium atom in a magnetic field using the Bethe- 
Salpeter equation without the equal-time Ansatz, survives the radiative corrections, as well 
as the effect of vacuum instability that occurs at the magnetic field value about b = 1.6 x 10^*. 
This indicates the existence of a maximum magnetic field in quantum electrodynamics. Note 



that contrary to the Dirac case 33|], where the critical magnetic field is determined by the 
large factor exp{l/aZ), for the Bethe-Salpeter case {Z = 1 for positronium) we got the 
factor exp(l/a;^/^). 

From Eq. (!67l) the relative correction to the ground state energy for moderate magnetic 
fields is 

Efl^^^L^V. (68) 

\Eo\ Stt \^ 2aZ J ^ ' 

For Z = 1 this correction, when extrapolated (though unrighteously) down to the value 
b = 0.27, is of the same order of magnitude (5.8 x 10~^) as the relativistic relative correction 



(2.8 X 10^^) calculated by Goldman and Chen 46| basing on the Dirac equation for this 
- largest in their analysis - value of b. The same situation retains, if the results of these 
authors are linearly extrapolated (using their two largest values of b for Z = 1) into the 
region of larger 6, 1 ^ 6 ^ 10^, wherein (1^51) is valid. Therefore, already for magnetic fields 
far from critical fields causing the free positron production instability the impact of vacuum 
polarization is at least no less important than relativism introduced by the use of the Dirac, 
instead of the Schrodinger equation. 
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VI. DISCUSSION 



In this paper we have shown that the electric field of a pointlike charge placed in a strong 
magnetic field {b = B / Bq ^ 1) may be significantly modified by the vacuum polarization, 
especially ii b > Sna^^ ~ 10'^. At present, it is commonly accepted that many compact 
astronomical objects identified with neutron stars are strongly magnetized. For soft gamma- 
ray repeaters and anomalous X-ray pulsars, for instance, the strength of the surface magnetic 
field is estimated as ~ 10^^-10^^ G 



47| . Several radio pulsars with similar surface magnetic 



fields have been recently discovered 



481. More strong magnetic fields {B ~ 10^*^ - 10^'^ G 



or even higher) are predicted to exist at the surface of cosmological gamma-ray bursters if 



they are rotation-powered neutron stars similar to radio pulsars |49j. The modification of 
the Coulomb law should affect the electric fields of an atomic nuclei and electrons placed 
in such a strong magnetic field. The electric field of a particle is one of its fundamental 
features. Therefore, at the surface of neutron stars with extremely strong magnetic fields 
many properties of matter (including individual atoms and molecules) and various physical 
processes (such as radiation of particles) where the electric field of particles is important 

50|) may be changed 



(for a review on physics of strongly magnetized neutron stars, see 
substantially by the present modification of the Coulomb law. One of such changes is 
discussed in Section V where we have come to negation of the standard result j32|, referred 
to in many speculations on behavior of matter on the surface of strongly magnetized neutron 
stars (e.g., [SOj and references therein), that the ground-state energy tends to negative 
infinity as the magnetic field unlimitedly grows. 

We hope that the modification of the Coulomb potential described in the present pa- 
per may lead to observational appearances in neutron stars with extremely strong mag- 
netic fields. As for the results relating to much larger magnetic fields, infinite in the limit, 
such as the QED string formation, these may be of fundamental importance as introduc- 
ing a nonempty magnetic-field-independent two-dimensional theory in virtue of dynamical 
dimensional reduction from 4-dimensional quantum electrodynamics. 
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Appendix I 

A. Asymptotic expansion around the singular point x = 

To consider the behavior of the potential near its pointlike source let us add to and 
subtract from (fT8|) the standard Coulomb potential ([5]) in the form 



' fr\?!!#)d.i.i^^ (69, 



so that 



ylo(x)=A^(x)-AAo(x), (70) 



where 

Note that the function Ay4o(x3, x_l) is an entire function of since the exponential in (fT5|) 
provides convergence of the integral ( 17Ti) for any complex value of this variable. Keeping 
quadratic terms in the power series expansion of Jo(^±a;_L) and exp(— 1/^3X3) in (|7T]l we 
obtain the first three terms of the asymptotic expansion of the potential (ITSi) near the origin 
0:3 = a;j_ = 



Ao(x) ~ 4^ - 2m(C - {2mx^YC^ - {2mx^YC\{)j , (72) 
where C, C± and C\\ are dimensionless positive constants depending on the external field: 



qm 7r2 \Am?J Jq (A;^ + A;|)(A;^ + A;| - /t2(0, A;^, /c^)) 
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"^11 SmTT^io Um'jJo {kl + kl){kl + kl-^,{0,klkl)r'' ^""^ 

Thanks to the exponential factor the integrals over k\ here are fast converging. The resulting 
functions decrease for large k^ as ~ 1/^1, so the remaining integrals over ks in fl73l) . fl7^ . 
(I75|) converge, bearing in mind that T is a bounded function. The inequality C± ^ C\\ 
implies the anisotropy. 

The values of the coefficients (1751) . (171|) . (1751) calculated for four values of the 
magnetic field h = W^, b = 10^, b = 10^ and b = 10^° are listed in the Table 



b 10^ 10^ 10^ 10^° 

C 2.21 9.08 31.37 32.70x10^ 
C± 75.9 2.58x10^8.38x10^ 8.49x10^° 
Cii 174.3 5.55x10^ 1.76x10^ 1.67x10" 



To find the asymptotic behavior of the constant C (!73|) as 6 — > oo, we may use first the 
representation fl24|) for the short-range part of the potential. The corresponding contribution 
Cs.T. into C is 

Cs.=M r(^- : " \ (76) 



By restricting the upper integration limit to the value y vr/a and substituting unity for the 



exponential the integral in fl76l) can be estimated as approximately equal to (1 — ^Ja/2Ti) = 
0.996. A computer calculation results in the value 0.9595. Correspondingly 

^ 0.9595^^. (77) 

The resulting values of C^,,, for b = 10^ 10^ 10^ 10^° are 3.27, 10.34, 32.70, 32.70x10^, 
correspondingly, to be compared with the exact values given in Table above. The co- 
incidence is the better the larger the field. It improves if the (negative) contribution 
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to C of the long-range part of the potential is added to to this row of numbers. (The 
absolute value of) the latter is a decreasing function of h that takes the limiting value 
Ci.r. = — (27r/gm)Ai r.(0, 0) |^^^ = —1.4152 according to Subsection B of Section IV. Note 
that if just the Yukawa law fl2^ is accepted for the potential we would deduce that for 



b = 10^, 10^ 10^ 10^° the values of C calculated following the Yukawa law are: 3.41, 10.78, 
34.01 and 34.08 x 10^. 

B. Modified Coulomb potential for less huge magnetic fields 

Consider the "moderate" values of the magnetic field in the interval 10'^ ^ 6 ^ 1, so 
that although b is large, but (a;6/27r) = 1.16 x lO^^b is still much less than unity. We shall 
present here the vacuum polarization correction to the Coulomb potential, which in this 
case is small. 

One may neglect k,2 in the denominator of (17T]) after the difference in it is completed and 
we obtain a magnetized vacuum analog of the Uehling-Serber potential Contrary to the 
Uehling-Serber potential that is of the order of a, its analog under consideration here is of 
the order of ab, i.e. much larger, given that 6^1. We shall be interested in x± = 0. Then, 
the /c_L-integral in AAo{x3,0) becomes 



Here Ei is the exponential integral, and we have used Eq. 3.353.3 from [41]. When integrating 
this over we may pass to the limit {2m?b/k1) = (2ei?/A;|) ^ oo in the integrand, since 
the remaining integral 



converges both at small and large integration variable (note the asymptotic behavior of 
( IT6|) ). Next we use the integral representation (fT6ll and the residue method to calculate 
the last integral. This leads to the additive vacuum polarization correction to the Coulomb 
potential in the form 




For the four values of the external field 




(78) 




(79) 




(80) 
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Setting X3 = in it we obtain for ([73] 



c.^. (SI, 



So, in the interval of magnetic fields indicated the constant C in the Laurent expansion f[72]) 
grows linearly with the field, in contrast to the square root growth (1771) characteristic of 
larger fields, as we saw in the previous subsection. The correction ( IHOi) was used in 35|] to 
find the energy correction to (l56ll for 1^6^ 1000. 



Appendix II 



In this Appendix we solve, for asymptotically large magnetic fields h ^ {2ti /a) ~ 10^, the 
eigenvalue problem inferred by the Schrodinger equation fj55|) with only the short-range part 
( !20l) of the modified Coulomb potential taken for Ao(x3, x_l = 0). The latter is approximated, 
in accord with ( 172|) with the quadratic terms omitted, = C\\ = 0, as eAs,r.(a;3, 0) = 



-Za ij^- 2mC 




for Lb < \xs\ < X3 
for |x3| > X3 = 



1 



2mC 

1 



(82) 



2mC ' 



where the external- field-dependent constant C ( !73|) is given by eq. (1771) . In the same way 
as in Sec. IIIB we may derive that the potential fl82l) becomes the 5-function in the b = 00 
limit, with the coefficient, however, different from the one in (1301) : 



V{xs) 



qe 
'2^ 



, Vb 
In — 

2C 



qe 
'2^ 



2aJ 



. 2C 
Vb 

TT 



S{xs) 

1/21 



5(X3) 



4.79|^5(X3). 

/TT 



(83) 



The difference in coefficients is owing to the fact that we kept only two terms in the expansion 
fl72]) . In equality ([HSD Eq- fl77|) was used. The square root asymptotic dependence fITTI) of C 
on the magnetic field is crucial for the formation of the 5-function limit of the potential. 

The approximation (1521) replaces the curves in Fig. [3] by continuous broken lines. The 
lowest energy state of the Schrodinger equation fl55|) is determined by imposing the boundary 



condition 



34| 



d^(a:3) 

dX3 



0. 



(84) 



X3=Lb 
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For the approximation (182|) to be meaningful it is necessary that 

xi > Lb. (85) 

With Eq. ( 1771) for C, this condition reduces to the evident inequality 1377r/2 ^ (0.9595)^ 
and is thus guaranteed. 

Introducing the so-called quantum defect 6 instead of the eigen-energy E according to 
the relation [remind that ab = {ma)~^ is the Bohr radius] 

^ - ^"^'"^ = -sJx- '''' 

and the new variable z = 2x3/ 6aB we obtain for ( |55l) two equations 

—\^ + —^{z)--^{z) = 0, for -^<z<z=—^ = — 87 

dz^ z ^ (5aB mC(5aB ^ ' 

and 

d2$(^) Zb'^C ^, ^ 1 w N ^ r 1 « 

— rV + ^0^-7^0^ =0, for z>z=— ^ = — . 88 

One should consider the couple of Eqs. ( 1871) and ( IHHil with the boundary condition 

d*(^) 



d^ 



2L£ 



(89) 



that follows from fl84j) . as an eigenvalue problem for determining the quantum defect 6 and 
hence the energ y (1861) . The general solution to the confluent hypergeometric differential 
equation (1571) is [51 1 



^ = ^W^z5,i(^)+M^5,i(^), ^<^<^, (90) 

where Wzgi{z) is the Whittaker function, decreasing at z 00, while the other, linear 
independent solution, growing at 2; — > 00, M^^i is expressed in terms of the confluent 
hypergeometric function $ as 

M^s,iiz)=e-h^l-Z6,2;z), (91) 

and A and B are constants. 

We shall seek for the solution of the eigenvalue problem (1871) . (!88l) . (189!) in the region 
[serving the asymptotically large magnetic fields considered here] 

6 » ^, (92) 
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so that z <^ 1. Therefore, only the small- distance behavior of the fundamental solutions to 



Eq. 

z 



3J, 



871) will be important. Referring to the asymptotic behavior of the solutions at small 
51| 

1 



exp — ^ 



+ ^[In z^^{\- Z6) - - tlj{2)] + 0{z^) In zj , 



(93) 



that retains the terms z^, z and 2; In 2; (here the logarithmic derivative ip of the Euler F- 
function F appears), Eq. fl^Ul) is matched continuously in the point z = z with the decreasing 
solution of the Schrodinger equation (|88ll 



z> z. 



(94) 



where 



i _ derivative over z 

4 a ' 



d^Jiz) 



Az 



-fs 



(95) 



if the coefficients A and B in (19011 are taken as 

1 „ „ 1 dWzsA^) 



A 



B = -h 



Wz5l{z) dz 



(96) 



Keeping the leading terms as 2; — > {Z5zhiz is neglected as compared to 1) we get from 



^Wzsi{z) 



2F(1 - Z5) 



+ 



In^ + 7 + ^(1 - Z5)-i){l) 
T{-Z5) ' 



(97) 



where 7 = — '?/'(l) is the Euler constant. With these values, the boundary condition ([89 
results in the following algebraic equation for the quantum defect 5 



= In 2; — In — ^ = — ln(2mLBC) . 
Zo oob 



(98) 



Solution to equation (l98l) is 



4^2 In^ (2mCLB) + 



4ZC\'^ 



a 



(99) 
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and the ground-state energy (l86l) is 

2^2 



E 



ma 



2 In^ {2mCLB) = -2Z^a'm\n^ ' ^ 



2C 



(100) 



It remains to make sure that the assumption (p2l) 1 = a/ C5 <^ 1 necessary for the use 
of the asymptotic form of the solutions (|93|) made while deriving equation fl98|) is indeed 
obeyed, once the quantum defect 5 is given by (!99l) . In other words, we must check the 
strong inequality 



In^f^V- 



1/2 ^ 

Za 



(101) 



By solving the quadratic equation with respect to C /Za it becomes 

1 + 



InM^ 



1/2 ^ 



2Za 



(102) 



For the fields so large that Eq. fl77|) holds for C, the In^-term becomes independent of 
the magnetic field 



In^ ( ^ 



TT 



ln%/ — = 7.213, 
^' 2a 



(103) 



and the inequality (11021) gives 



6 > 1 + W 1 + n , , 

0.9 1 V V 2a 




MaZHQ.Q = 3Z^ 



(104) 



This condition on the values of the magnetic field is less restrictive than the condition 
of validity of Eq. ( iTTl) . Thus the inequality ( !92l) is aposteriory verified. Note, that also 
the inequality ^5 <^ 1 is satisfied for the same fields, justifying the disregard of ZS'z In^ 
made when writing ( 1971) . We conclude that for asymptotically strong magnetic fields the 
derivation that has led to Eq. (llOOp is justified, and the ground-state energy acquires the 
magnetic-field-independent limiting value 



E = -2Z^a^m\n\l — = -7.686 x lO'^mZ^ = -389.3 eV x Z' 



TT 




2a 



(105) 



Note added in proof. Most recently a work 52] appeared, where the modified Coulomb 



potential in a strong magnetic field calculated in |35|, l36| and in the present paper is also 



considered. In particular, in that work the long-range asymptotic behavior |x| given as 
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Eq. (138|) is supplemented by two next-to-leading terms of the order of |x| ^ and |x| 



Besides, the authors of 



52j found a small (~ a/vr) anisotropic amendment to the simplified, 



Yukawa-like, form of the exact scaling equation 
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